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ABSTRACT 

In  the  present  work ,  structural  analysis  of  stiffened  panel  for  aircraft  fuselage  is  performed  for  pure  composites 
and  composite  plates  with  isogrid  and  orthogrid  using  FEM.  In  the  bulking  analysis  of  the  panel ,  the  static  analysis  is 
carried  out  under  the  specified  loading  and  boundary  conditions  and  results  are  extended  to  the  buckling  analysis.  In  the 
buckling  analysis  of  stiffened  panel ,  the  mode  shapes  are  extracted  and  frequencies  are  tabulated.  The  material  selected 
for  the  analysis  is  carbon  fiber ,  carbon  fabric  epoxy  composite  by  precision  contact  lay-up.  The  fiber  orientation  for  the 
isogrid  structure  is  60°/30°  and  for  the  ortho  grid ,  the  fiber  orientation  is  0°/90°.  The  analysis  shows  that  the  orthogrid 
structure  is  better  than  the  other  panels.  The  inplane  load  bearing  capacity  of  the  orthogrid  structured  panel  is  more  than 
the  isogrid  panel  structure.  The  weight  of  the  orthogrid  structured  panel  is  less  than  the  isogrid  stiffened  panel. 
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1.  INTRODUCTION 

Fuselage  is  a  long  hollow  tube  that  carries  payload  and  the  equipment  required  for  the  sustainable  flight. 
The  aircraft  fuselage  is  made  up  of  aluminum  rings  and  stiffeners  to  support  a  thin  aluminum  skin  to  reduce 
weight.  An  aluminum  fuselage  is  typically  constructed  by  the  semi-monocoque  method,  which  consists  of 
longitudinal  stringers,  ring  frames,  and  end  rings  as  shown  in  Figure  1.1.  A  series  of  frames  in  the  shape  of  the 
fuselage  cross  section  are  held  in  position  on  a  rigid  fixture.  Aluminum  ring  frames  are  equally  spaced  over  the 
length  of  the  cylinder.  The  ring  frames  are  joined  with  light-weight  longitudinal  elements  called  stringers.  The  ring 
frame  and  stringers  are  covered  with  a  skin  of  sheet  aluminum,  which  is  attached  by  riveting  or  by  bonding  with 
special  adhesives. 


Figure  1.1:  Fuselage  Showing  Frames,  Stringers  and  Skin 


1.1.  Fuselage  Panel 

This  is  a  simple  structural  analysis  of  a  composite  fuselage  panel  of  a  model  airplane  made  of  stringers, 
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i.e.  a  stringer  riveted  to  a  skin  panel  [1].  The  fuselage  section  is  of  uniform  configuration  along  its  length  and  its  cross- 
sectional  area  is  defined  to  be  a  plate  as  shown.  Figure  1  shows  the  whole  panel  and  gives  an  idea  of  its  location  in  a  real 
aircraft.  The  considered  stringers  have  a  trapezoidal  profile  -  they  are  also  called  Omega  stringers.  It  is  held  fixed  to  the 
body  of  the  airplane  on  one  end  and  pinned  at  the  other  end.  The  panel  and  the  associated  super- stringers  are  made  of 
composite  layers,  each  of  them  being  defined  by  an  orientation  and  a  thickness 


Figure  1.2:  Fuselage  Panel  and  Associated  Model 


Figure  1.3:  Omega  Super-Stringer  (Stringer  Profile  and  Adjacent  Skin  Panels) 
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Figure  1.4:  Laminated  Composite  Plate 


While  the  fiber  angle  is  restricted  to  take  eleven  discrete  values(90°,30°,-30°,  30,  -30,  -30,  30,  -30,  30and90°  ),  the 
thickness  may  vary  continuously  between  some  predefined  bounds,  which  were  set  to  0.5  and  15  mm  in  this  case.  The  goal 
of  our  problem  is  to  obtain  Buckling  behavior  of  given  fuselage  sections,  where  the  first  buckling  loads  is  evaluated  and 
corresponding  mode  shapes  will  be  presented,  and  to  find  which  comparing  the  result  of  each  plate  for  its  efficiency  using 
finite  element  analysis  method. 


1.3.  Conventional  Fuselage  Loads 

The  fuselage  structure  design  should  sustain  the  cyclic  pressure  loads  due  to  cabin  pressurization  and  De- 
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pressurization,  Bending,  Shear,  Torsional  and  fatigue  loads.  The  fuselage  should  be  designed  for  the  strength,  rigidity  and 
fail  life. 

The  direct  stress  distribution  due  to  bending  can  be  calculated 

Af  (/  v  — |  I  0  A t  (f  x  -  I  y) 

fT  =  - ; - 1 - ; - : - 


Where, 

MXj  My  =  Bending  moment  about  x  and  y  axis 

Ixx,Iyy,  Ixy=  Area  moment  of  inertia. 

The  shear  flow  distribution  due  shear  can  be  obtained  as 

( $  f  -Si  ( S  i  -Si  i  ■ 

^  Bf  y  -  - V2  T  +  q3 


Where, 

Qs=  Shear  flow  distribution 
Sx,Sy=  Shear  loads  in  x  and  y  direction 
Ixx,Iyy,Ixy=  Moment  of  areas 
Xr,Yr  =  Center  of  twist  at  point  x,  y 
Br=  Cross-sectional  area 

The  shear  flow  due  to  Torsion  for  single-cell  closed-section  beam  can  be  obtained 

r 

where 

q=  Shear  flow  distribution  due  to  applied  torque 
T=  Applied  torque 
A=  Area 

2.1.  DESIGN  OF  STIFFENED  FUSELAGE  PANEL 
Isentropic  Material 

The  panel  of  the  fuselage  section  is  made  from  the  isotropic  aluminum  plate.  The  dimensions  of  the  panel  taken 
for  the  analysis  is  shown  in  the  figure2.1. 
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Figure2.1:  Isotropic  Plate 


Table  1.1:  Thickness  of  the  Stiffened  Panel 


Material  Type 

Skin  Thickness 

Stiffner 

No  of 

Thickness 

(mm) 

Height 

Layers 

Isotropic(aluminum) 

7 

- 

- 

6.90  mm 

Pure  composite 

5 

15 

10 

13 

Composite  with  isogrid 

5 

15 

10 

20 

Composite  with  orthogrid 

5 

15 

10 

20 

Thickness  Calculation 

The  plate  is  fixed  at  one  end  and  hinged  at  other  end  as  it  is  a  long  plate  made  up  of  isotropic  material. 
Given  Data:  L=  800  mm 
B  =  200  mm 
Operating  load  =  2  tons 
Factor  of  safety  =1.33 

Buckling  factor  =design  load  /  Operating  load 


Figure2.2:  Plate  with  Loads  &  Constraints 

Design  load  =  buckling  factor  x  operating  load 
=  5  x2000x9.81  x  1.33 
=  130470  N 

Design  load  (Pe)  =  130470  N 
The  buckling  load  is  obtained 


Pe  =tt2EI  /Le2 

130470  =(ti2  x  210  x  103  x  150  x  T3)  /  (8002/(V2)2  xl2) 
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t  =  6.90  mm  (approx  7  mm  ) 

The  thickness  for  a  isotropic  plate  is  t  =  7  mm 
Composite  Plate 

The  fuselage  panel  is  made  up  composite  material 


Figure2.3:  Composite  Plate 

Dimensions  of  Composite  Plate 

Thickness  of  plate  =t  =  13  mm 

Thickness  Calculation 

For  a  long  plate  with  all  sides  clamped  and  loaded  in  the  x-  direction  only,  the  buckling  load  per  unit  width  can  be 
estimated  as 

NxCr=(  n2/t>2  )  (4.6  (  DllD22)1/2  +  2.67  D12  +  5.33  D66) 

Substituting  the  value  of  constants 

26600  =  (  n2/1502  xlOOO)  (4.6  (  6.0605T3  x  1.322T3)!/2  +  (2.67  xl.87)  +  (5.33x2.056)T3 
t  =  12.77  mm 

The  thickness  of  pure  composite  plate  =  t  =  13mm 
Stiffened  Panel  with  Isogrid 


Figure2.4:  Composite  Plate  with  Isogrid  Stiffner 

•  Dimension  of  composite  plate  with  isogrid  stiffner 

•  Length  =  800  mm 

•  Breath=  200  mm 

•  Skin  thickness^  5  mm 
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•  Stiffner  thickness=  15  mm 

•  Angle  between  stiffner s  =  60°  and  30° 

Total  numbers  of  layers  are  eleven  out  of  which  first  ten  layers  are  of  skin  and  remaining  one  is  of  stiffner  having 
thicknessof  15  mm 

Each  layer  from  one  to  ten  is  of  thickness  0.5  mm 

Total  thickness  of  plate  =(skin  thickness  +  stiffner’ s  thickness) 

=  (0.5  x  10)  +  15) 

T  =  20  mm 

Stiffened  Panel  with  Orthogrid 


A  it 


Figure2.5:  Composite  Plate  with  Orthogrid  Stiffner 


Dimensions  of  composite  plate  with  orthogrid  stiffner 

Length^  800  mm 

Breath  =  200  mm 

Skin  thickness^  5  mm 

Stiffner  height  =15  mm 

Angle  between  stiffner’ s  =  0  and  90  degree 

Total  numbers  of  layers  are  eleven  out  of  which  first  ten  layers  are  of  skin  and  remaining  one  is  of  stiffener  having 
height  of  15  mm  Each  layers  from  one  to  ten  is  of  thickness  0.5  mm. 

Total  thickness  of  plate  =  (skin  thickness  +  stiffner  height) 

t=  (0.5  x  10)  +  15) 

t  =  20  mm. 

3.1:  Buckling  Analysis  of  the  Fuselage  Panel 

The  stiffness  matrix  of  an  element  [K]  is  calculated  from  the  equation, 

K  =  }v  [B]T  [C]  [B]  dv 
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Where,  [B]  is  the  strain-displacement  matrix. 

[C]  is  constitutive  matrix. 

Formulation  of  Geometric  Stiffness  Matrix  [Kg] 

Geometric  stiffness  matrix  can  be  written  as 
[Kg]  =  Jv  d[Bnl]T{c}dv 

Where  [Bnl]  is  the  nonlinear  strain-displacement  matrix,  {a}  is  the  compressive  stress  matrix 

The  global  stiffness  matrix  of  the  stiffened  plate  finite  element  model  is  obtained  by  computing  the  element 
stiffness  matrix  of  each  element  and  assembling  them  by  posting  them  in  appropriate  global  locations  determined  by  node 
numbering  and  connectivity. 

Linear  Buckling  Analysis 

During  buckling,  the  total  stiffness  matrix  becomes  singular  or  the  determinant  of  the  total  stiffness  matrix 

vanishes. 

The  eigen  value  problem  of  instability  is  therefore  formulated  as  ([Ko]  +[KG])  {8}  =  0  (1.2) 

([Ko]+5b  [Kg])  {8}=  0(1.3) 

The  buckling  pressure  is  evaluated  for  the  condition  [Ko]+  Ab[Kg]  =  0  (1.4) 

Where,  b  is  the  non-dimensional  buckling  pressure. 

The  solution  procedure  adopted  for  obtaining  the  linear  buckling  pressure  is  the  determinant  search  procedure. 

A  value  of  8  is  assumed  and  the  determinant  of  the  matrix  is  calculated. 

The  process  is  repeated  by  Element  stiffness  matrix  [K]: 

The  stiffness  matrix  of  an  element  [K]  calculated  from  the  equation, 

K  =  jv  [B]T  [C]  [B]  dv  Where  [B]  is  the  strain-displacement  matrix.  [C]  is  constitutive  matrix. 

Formulation  of  Geometric  Stiffness  Matrix  [Kg] 

•  Geometric  stiffness  matrix  can  be  derived  from  the  expression  [16]  given  below. 

•  [Kg]  =  jv  d[Bnl]T{a}dv  Where  [Bnl]  is  the  nonlinear  strain-displacement  matrix,  {a}  is  the  compressive  stress 
matrix 

•  Assembly  of  global  stiffness  matrices  and  application  of  boundary  condition 

•  The  global  stiffness  matrix  of  the  stiffened  plate  finite  element  model  is  obtained  by  computing  the  element 
stiffness  matrix  of  each  shell  element  and  assembling  them  by  posting  them  in  appropriate  global  locations 
determined  by  node  numbering  and  connectivity. 

•  Linear  buckling  analysis  During  buckling  the  total  stiffness  matrix  becomes  singular 

•  The  eigen  value  problem  of  instability  is  therefore  formulated  as  ([Ko]  +[KG])  {8}  =  0(2) 
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•  ([Ko]+  Xb  [Kg])  {8}  =  0(3) 

•  The  buckling  pressure  is  evaluated  for  the  condition  [Ko]+  Xb[Kg]  =  0(4) 

The  finite  element  method  has  become  a  powerful  tool  for  the  numerical  solution  of  a  wide  range  of  engineering 
problems.  An  assembly  process,  duly  considering  the  loading  solution  of  these  equations  gives  us  the  approximate 
behavior  of  the  continuum.  ANSYS,  Inc,  is  an  engineering  simulation  software  provider  founded  by  software 
engineer  John  Swanson.  It  develops  general-purpose  finite  element  analysis  and  dynamics  software.  While  ANSYS  has 
developed  a  range  of  computer-aided  engineering  (CAE)  products,  it  is  perhaps  best  known  for  its  ANSYS  Mechanical  and 
ANSYS  Multi  physics  products. 

4.  RESULTS  AND  DISCUSSIONS 


The  buckling  analysis  of  the  Iso-grid  and  Ortho-grid  stiffened  panel  structure  has  been  performed  in  this  work. 
Based  on  the  results  obtained  from  the  FE  ANSYS  analysis,  the  Ortho-grid  stiffened  structure  is  better  than  the  Iso-grid  in 
terms  of  the  load  bearing  capability  and  weight.  The  deformation  for  different  materials  under  given  loading  conditions  are 
tabulated.  The  deformations  for  the  isogrid  and  pure  composite  are  more  than  the  other  materials.  Similarly,  the  frequency 
for  different  materials  on  the  same  loading  and  boundary  conditions  are  shown  in  the  in  the  table.  The  frequency  for 
composite  with  isogride  is  lower  than  the  frequency  of  other  materials.  The  Configuration  of  fuselage  section  is  shown  in 
the  Figure.4.1  and  the  Figure.4.2  depicts  the  panel  taken  for  the  analysis. 


Figure.4.1:  Orthogrid  Fuselage  Section  Figure.4.2:  Orthogrid  Stiffened  Panel 


Figure  4.3:  Ply  Sequence  at  Skin  Figure  4.4:  Load  &  Boundary  Conditions 
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Figure  4.5:  Load  &  Boundary  Conditions  Figure  4.6:  Buckling  Mode  Shape 


Table4.1  frequency  and  Deformation 


MATERIAL  TYPE 

Load(kN) 

Deformation  (mm) 

Frequency  max  (Hz) 

Isotropic(aluminium) 

2660 

0.01403 

6.849 

Pure  composite 

2660 

0.024961 

4.34 

Composite  with  isogrid 

2660 

0.021614 

3.905 

Composite  with  orthogrid 

2660 

0.017179 

4.671 

5.  CONCLUSIONS 


The  stiffened  panel  for  the  fuselage  has  been  designed  and  analyzed  for  Buckling  using  Ansys  software.  The  level 
of  deformations  and  frequencies  observed  are  close  to  the  theoretical  values.  However,  the  results  of  the  laminate  analysis 
and  finite  element  analysis  shows  that  they  are  within  10%.  The  maximum  deformation  of  the  composite  panel  under  given 
loading  and  boundary  conditions  is  0.017179mm,  which  is  very  less  compared  to  other  materials.  The  experimental  results 
of  stiffened  panels  for  buckling  response  under  compression  loading  were  compared  with  the  numerical  results  obtained 
using  linear  and  nonlinear  finite  element  analysis. 
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